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■ Abstract. We describe an optomechanical system in which the mean phonon number 

of a single mechanical mode conditionally displaces the amplitude of the optical field. 
' Using homodyne detection of the output field we establish the conditions under which 

phonon number quantum jumps can be inferred from the measurement record: both 
the cavity damping rate and the measurement rate of the phonon number must be much 
' greater than the thermalization rate of the mechanical mode. We present simulations of 

^ , the conditional dynamics of the measured system using the stochastic master equation. 

In the good- measurement limit, the conditional evolution of the mean phonon number 
shows quantum jumps as phonons enter and exit the mechanical resonator via the 
bath. 
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1. Introduction 



Observation of energy quantization in a macroscopic mechanical oscillator mode is of 
fundamental interest in exploring the quantum-to-classical transition. Electromechani- 
cal systems, where the relevant mechanical mode interacts with devices such as super- 
X : conducting quantum interference devices, single-electron transistors, cooper-pair boxes, 

, or even other mechanical modes, provide a possible route toward achieving this, and 
various scenarios have been explored theoretically [I]-[8]. Experimental realization, 
however, remains a challenge. A second platform where observation of macroscopic 
mechanical mode energy quantization may become possible is optomechanics [9l [TO] . 
In an optomechanical system, a collective vibrational mode of a mesoscopic system, 
such as a doubly clamped beam |Tl], a ring [12], a suspended membrane ("membrane- 
in-t he-middle" ) [131 111] 5 or even a cloud of ultracold atoms [I5]-[T7], is parametrically 
coupled to an electromagnetic cavity field mode. While also challenging, there has been 
much interest in the prospects of experimental observation of macroscopic mechanical 
mode energy quantization in an optomechanical system due to recent developments. In 
particular, direct coupling of the optical cavity mode to the square of the mechanical 
displacement coordinate has been experimentally demonstrated in the " membrane-in- 
the- middle" system [131 E] its ultracold atoms analogue [E]. Under the rotating 
wave approximation, such a coupling allows for a quantum nondemohtion (QND) [18] 
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measurement of the mechanical mode phonon number. Observation of mechanical en- 
ergy quantization in such a system can then be achieved either indirectly, through 
phonon shot noise measurement, or directly, by detecting quantum jumps in the phonon 
number via continuous measurement of the cavity field. A. A. Clerk et al. have analyzed 
the former case in [19]. In this paper we present an analysis of the latter. In partic- 
ular, we show the nature of the displacement of the cavity field in the strongly driven 
regime in response to changes in the phonon number, and we establish constraints on 
the parameter space which must be satisfied in order for quantum jumps of the phonon 
number to be observed. 

There are two ways we can describe a measurement in quantum mechanics. Firstly 
we can simply give the solution to the master equation describing the interaction 
between the measured system and the apparatus and any decay channels that may be 
present. This is called the unconditional evolution. Secondly, in the case of a continuous 
measurement we wish to describe the conditional state of the measured system given a 
particular classical stochastic record for the measurement. For the case considered here 
this will require a phase-sensitive measurement of the cavity field, such as homodyne 
detection. The measurement record is then a stochastic homodyne current and the 
evolution of the mechanical resonator is conditioned on a particular realization of this 
current [18] . Of course the statistics of the homodyne current itself is partly determined 
by the quantum state of the mechanical resonator. The goal is then to reconstruct 
the quantum trajectory of the mechanical mode by suitably filtering the measurement 
record, and to observe phonon absorptions and emissions of the mechanical resonator 
mode as retroactive quantum jumps [20] in the trajectory. If one averages the conditional 
state of the system over all possible measurement records, one obtains the unconditional 
state. 

2. The model 

We assume the case of a nearly ideal single-port cavity such that the constraint 
established by Miao et al. in [21] is satisfied. Under a single-mode description of 
the cavity the generic Hamiltonian for the optomechanical system with a drive on the 
cavity is 



where a, a^" refer to the optical mode and 6, 6^ refer to the mechanical mode, uJc{x) 
is the cavity mode frequency dependent on the mechanical resonator position x, z/ is 
the frequency of the mechanical mode, u and e refer respectively to the frequency and 
amplitude of the coherent drive field on the cavity, Xq is the equilibrium position of the 
mechanical resonator, and 



(1) 




(2) 
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is the displacement coordinate of the mechanical mode. If the mechanical resonator is 
positioned such that xq is an extremum of ujc{x), we approximate 

H = hucixo)a^a + hub^b + h{e*ae''^' + ea^e"*"'*) + ^^0^0(6 + b^f, (3) 



where 



2i^m dx"^ 



(4) 

x=xa 

We now assume that the optical cavity remains near the steady state it would 
have reached if G = 0. In an interaction picture at the driving frequency, this is a 
coherent state with amplitude ao = — ' where 5 = Uc — oj is the detuning between 
the cavity driving field and the cavity frequency, and k is the cavity damping rate. 
With no loss of generality we will assume that the phase of e may be chosen so that 
ao is real. We now expand the interaction around this steady state value by making 
the canonical transformation a = a — a^. After the rotating wave approximation, the 
effective Hamiltonian in the interaction picture may then be written as 

Hi = ^x{a + -a^)b\ (5) 

where 

X = 2Gao, (6) 

and we have taken 5 = for simplicity. If we now include the damping of the cavity and 
the mechanical resonator, the total system state is described by the master equation 

^ = ~[Hj, p] + KV[a]p + 7(iV + l)V[b]p + -fNV[b% (7) 

where T>[c\p = cpc^ — c^cp/2 — pc^c/2, N is the mean thermal occupation of the 
mechanical resonator bath at frequency z/, 7 is the mechanical mode damping rate, 
and where we have dropped the bar on the d for simplicity of notation. 

The interaction Hamiltonian in equation indicates that the phonon number 
operator for the mechanical resonator is not changed by the interaction with the optical 
field and, in the absence of dissipation, is a constant of the motion. The model thus 
describes a QND measurement of the mechanical resonator's phonon number with the 
cavity field forming the first stage of the measurement apparatus. A similar model was 
considered in Walls et al. |22] for both components realized as optical field modes, but 
only mode-a was damped. 



3. Unconditional dynamics 

3.1. Without mechanical damping 

As pointed out in [2], for direct phonon number detection the measurement apparatus 
must have a dynamical time scale much shorter than that of the phonon number of 
the relevant mechanical mode. Therefore we begin by considering the dynamics on a 
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time scale over which the mechanical damping can be neglected. Such a problem is 
represented by the master equation above with 7 = 0. The solution to this problem was 
given by Walls et al. in |22] : 

,2 



p{t) = ^ Pnm{a,a')exp 



k2 



(n 



m] 



[I - Kt/2 - e-"*/2) 



xexp 



[n 



m] 



—i—a * — i—a I (1 

K K 



-Kt/2\ 



(l«n(t))«(t)|)(a) 



{a'n{t)\an{t)) 

where Pnm{c(,C(') is the initial probability distribution such that 

(l«)(«'l)(a) 



p(0) = J2 Pnm{a,a'){\n){m\)^b)® 



{a'\a) 



(8) 



(9) 



and and a'm{t) are complex amplitudes corresponding to coherent states of the 

cavity with a and a' as the initial amplitudes: 



Xn 

an{t) = -t — (1 

K 



ae 



1 _ e-«*/2) + 



(10) 
(11) 



K 



The solution in equation ([8]) illustrates some important features of quantum 
measurement. Firstly, for sufficiently strong coupling, the density operator rapidly 
becomes diagonal in the number basis for the resonator, which is the measured quantity. 
Secondly, the optical system is driven towards coherent states which are the pointer 
basis states [23j for the damped cavity field. The rate of diagonalization is proprotional 
to x^ I n. The resulting mixture is a classical correlation between Fock states of the 
resonator and corresponding pointer basis states for the field, and we find in the steady- 
state (a)(t) = —i-Ub, where nt, = (b^b). 



3.2. With mechanical damping 

Unravelling the full master equation with nonzero 7 is intractable, but since we are 
required to be in the adiabatic limit we can build our intuition of the full dynamics 
from the previous section. The nature of the cavity-environment and cavity-mechanics 
interaction is the same for both zero and nonzero 7. We therefore expect that the 
classical correlation between mechanical Fock states and cavity coherent states is 
preserved. It is only the initial distribution of the mechanical Fock states that will 
evolve and the cavity state will adiabatically follow so that (a)(t) ~ —i^ni,{t). We 
can verify this by solving the moment equations that follow directly from the master 
equation ([71) (which correspond to moments of the unconditional state). The moment 
equations are 
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The solutions are 

nfe(t) = n,(0)e-^* + N{1 - e"^*), (14) 

(15) 



(a)(t) = (a)(0)e-''*/2-.| 



In the adiabatic hmit we recover {a)(t) ^ ~i^nb{t), as expected. It is then apparent 
that we can extract information about the average phonon number by monitoring the 
quadrature phase amplitude to the steady state amplitude ao (which has been chosen 
as real). One expects that, prior to the optomechanical coupling turning on, the 
mechanical resonator will be in thermal equilibrium with its environment, in which 
case nb(0) = N. For long times, the change in the cavity amplitude from the steady 
background amplitude, ao, Aa = —ixN/n. We thus regard the ratio x/^ as the gain 
of the measurement. 

As demonstrated by equation f lT4|) . phonon number quantum jumps are not present 
in the unobserved system even in the quantum limit. Jump-like behaviour arises only 
in the presence of a continuous measurement, which is the subject of the next section. 



4. Conditional dynamics 

We have shown that information about the phonon number is reflected in the phase 
quadrature amplitude of the cavity field. Under continuous homodyne measurement 
of this quadrature, the system is governed by the following stochastic master equation 
(SME): 

dp = -l[Hj, p]dt + 7(iV + l)V[b]pdt + -fNVlb^pdt + KV[a]pdt 
n 

+^dWH[ae-'^p. (16) 

Here, T>[c\p is as defined for equation ([7]), 'H[c]p = cp + pc^— Tr(cp + pc''") is 
the measurement superoperator, rj is the detector efficiency, and dW is the Wiener 
increment. The homodyne measurement signal will be a photocurrent proportional to 
the cavity phase quadrature amplitude with stochastic noise due to the local oscillator 
and intrinsic quantum noise of the cavity field: 

ih{t) = r]K{ae''^ + a^e'^) + ^i{t), (17) 

where ^(t) = dW/dt. As previously mentioned, we require the cavity to adiabatically 
follow the mechanical number state. In this limit, we may adiabatically eliminate the 
cavity field from the full SME in order to reduce computational overhead. Solving the 
quantum Langevin equation for the steady state of the cavity field gives a = —i^b^b. 
The resulting SME for the density matrix of the mechanics is 

dph = 7(iV + l)V[b]pbdt + -fNV[b^]pbdt + TV[b^b]pbdt + ^/^H[b^be'''']pbdW, (18) 

2 

where F = From this it is straightforward to further simplify to the diagonal elements 
of pfe: 

dp„ = 7iV - (n + 1 )p„] rft + 7 (iV + 1 ) [(n + 1 )p„+i - J c/t - 2 v^(n - (ri) )M1^. ( 1^ 
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In the adiabatic limit the photocurrent becomes 

ih{t) = -2rjx{b^b), + (20) 

From hereon we set 77 = 1. We note that equation fll9p is of the same functional form 
as the one dealt with in Section V. of [2], and we refer the reader to that reference for 
complementary analysis and discussion. 

Under continuous measurement the evolution of the system is the result of a 
competition between two dynamical processes. The correlation between cavity coherent 
states and mechanical phonon number established in the previous section indicates that 
measurement of the cavity field will tend to collapse the variance of the mechanical 
state in the number state basis. The thermal bath on the other hand will have the 
opposite effect. For jump-like behaviour of the phonon number to arise, the measurement 
rate must dominate the thermalization rate. This amounts to the adiabatic and fast- 
measurement conditions, which we now define. 

In the absence of measurement, the thermalization rate of a mechanical Fock state 
\n) can be shown to be -i[N{n + 1) + (A^ + l)njil [2J. To be in the adiabatic limit, we 
therefore require the adiabatic condition: k ^ 'y[N{njnax + I) + {N + l)nmax], where 
iT-max is the largest phonon number state one wishes to resolve quantum jumps from. 

The fast-measurement condition arises due to the fact that even in the adiabatic 
limit the measurement will extract information on the phonon number at a finite rate. 
The rate of number state collapse for a system described by equation ( fT9l) is of the order 
r [21 121]. The fast-measurement condition is therefore F ^ 'j[N{nmax + '^) + {N+l)nmax]- 
This condition, along with the adiabatic condition, defines the "good-measurement 
limit," which we explore below. We note that this limit is valid for arbitrary values 
of the ratio x/^; the strong-coupling regime (x/^t ^ 1) is not required. 



4.1. Good-measurement limit 

In the good-measurement limit the cavity adiabatically follows the number state of the 
mechanical mode, and the measurement collapses the state of the mechanical mode in 
the number state basis very rapidly. Therefore we expect that the mechanical state 
will be a nearly pure number state most of the time with stochastic excitations 

(n — )• n-t- 1) or decays (n — )■ n — 1) due to the thermal bath. From the analysis of section 
III., the corresponding cavity field will have 

(a) = an{t), (21) 
with a quantum jump in the phonon number, 

\n){n\ — >\n±l){n±l\, (22) 
at time tj reflected as 

anit) ^ a.±i(t) = - e--(*-*^)/2) + a„(t,)e--(*-*^)/^. (23) 

I This expression is an upper bound on the thermahzation rate; under continuous measurement the 
thermahzation rate is decreased due to the measurement-induced coUapse. 
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The adiabatic condition ensures that the cavity amphtude will reach its steady-state 
value, —i^n, between consecutive jump times tj, and the phase quadrature amphtude 
will therefore trace out a step- like trajectory. 

To verify this intuition, we numerically integrate equation ( IT6|) with N = 0.5, 
x/k = 1.5, K = 1007iV, and x^/f^ = so that we are in the good-measurement 

limit for low phonon numbers. We start with the mechanics in the ground state. The 
result, figure [H shows well-resolved quantum jumps in the phonon number replicated 
by the cavity phase quadrature amplitude with a gain of x/'^ = The small 

variance of the phonon number distribution indicates a high purity of the mechanical 
state. The quantum trajectory formalism that we use assumes the ideal limit of 
an infinite amplitude local oscillator, thereby giving the photocurrent white noise 
with infinite amplitude in the limit that dt ^ 0. Numerically, the finite time step 
limits the noise amplitude and bandwidth. This is qualitatively consistent with the 
experimental situation of a finite amplitude local oscillator and a photodetector with a 
finite bandwidth. Still, the raw measurment signal, equation [T71 is likely to be dominated 
by noise from the local oscillator. However, we show below in section 4.4 that a simple 
sliding time-average of the raw signal is sufficient to reveal the quantum jumps. We plot 
such a low-pass filtered version of the homodyne photocurrent in figure [T]^d). 
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Figure 1. Evolution of: (a) phonon number, (b) phonon number variance, and (c) 
cavity field phase quadrature in the good-measurement limit for low phonon numbers 
{n ^ 1). Quantum jumps in the phonon number are replicated by the cavity phase 
quadrature amplitude, (d) The filtered homodyne current gives a noisy version of the 
cavity phase quadrature trajectory. 
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4-2. Verification of adiabatic condition 

The SME of equation ( fT6l) is valid outside of the adiabatic hmit and we can use it to 
verify the necessity of the adiabatic condition (k ^ 'y[N{nmax + 1) + (iV + l)nmax])- 
We set x^/i^ = 10^7iV and N = 0.5 so that the fast-measurement condition is satisfied 
for low phonon numbers, and perform simulations with k = 'jN, lO'-fN, and 10^7iV. 
The mechanical mode is initialized in the ground state. Trajectories of (6^6) for these 
simulations are shown in figure O We see that quantum jumps only arise when the 
adiabatic condition is satisfied, and thereafter are increasingly well-resolved for larger 
values of k. This simulation of figure [2]^a) has a shorter run time than figures [21(b) 
and (c) due to computational limitations, but the expected steady-state is reached: the 
measurement is too weak to induce number state collapse, and the mechanical mode 
thermalizes so that fib = N. 




Figure 2. Verification of adiabatic condition. Tlie trajectories slrow the evolution 
of tlie average phonon number with parameters N = 0.5, x^/^ = lO'^^N, and: (a) 
K — 77V, (b) K — lO'yN, (c) k — lO'^'^N. Jump-like behaviour arises only when 
K 3> j[N{n + 1) + {N + l)n\, where n is the phonon number. 



4-3. Verification of fast-measurement condition 

To verify the fast-measurement condition {x^/k ^ l[N{nmax + 1) + (iV + l)nmax]) 
we numerically integrate equation f lT9|) with k, = 10^7iV for three cases: x^/^^ = 
X^/k = 107iV, and x^/'^ = 10^7iV. We start with the mechanics in the ground state, 
and the bath temperature is set at iV = 0.5 so that the adiabatic condition is strongly 
satisfied for phonon numbers close to zero. The first case, figure |3]^ a), does not satisfy 
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the fast-measurement condition and therefore does not resolve quantum jumps in the 
phonon number. The second case, figure [3](b), is on the border of the fast-measurement 
regime for n ~ 1 and shows some jump-hke behaviour in the phonon number. The 
third case, figure |3]^c), strongly satisfies the fast-measurement condition for low phonon 
numbers and shows well-resolved quantum jumps in spite of being deeply within the 
weak coupling regime with x/k = 10~^. 




Figure 3. Verification of fast-measurement condition. The trajectories show the 
evolution of the average phonon number with parameters k — 10'*7-/V and: (a) 
X^/k = 7iV, (b) x^/k = 107iV, (c) / n — IQ'^^N. Jump-like behaviour occurs 
only when / n 3> ^[N(n -|- 1) -|- {N + 1)^1], where n is the phonon number. 

4.4- Processing the measurement record 

We must consider that although the adiabatic and fast-measurement conditions are 
sufficient for phonon number quantum jumps to arise in the phonon number trajectory, 
that trajectory must be inferred through the measurement record, which contains 
additional noise. Integration of the measurement record over a finite time interval 
can reduce the noise, but the time interval must be much less than the typical lifetime, 
'Tn^ao: ~ {liNirimax + I) + {N + l)nmax])~^ , of Fock State \nmax) m order for the jump- 
like behaviour between Fock states {|n) : n < Umax} to be resolved. Here we show that 
simple averaging of the measurement record for a finite time interval 5t ^ Tmax always 
allows one to infer the phonon number trajectory with resolution sufficient to see the 
quantum jumps. 
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We define a random variable x which is the integral of the homodyne current over 
this interval of time: 

X = ih{t)5t = -2x{b^b),6t + y^dW{t). (24) 

This is equivalent to a generalized measurement of the number described by the 
conditional probability distribution 

{x + 2xn5tf 



2K6t 



P(x\n) = (2-KK6t) ^^^exp 
The probablity distribution to get the result x is given by 

oo 

P{x) = Y,Pn{t)P 



[x\n] 



(25) 



(26) 



n=0 



where Pn{t) is the phonon number distribution at the start of the time interval 6t. One 
can then easily verify that 

X = -2x{n)6t, (27) 
var{x) = ^x^St^ var{n) + K^t, (28) 

so that the added noise is diffusive as expected from equation (1241) . The a posteriori 
state of the mechanical resonator given a particular result x is given by Bayes' rule as 

Pn{t)P{n\x) 



P^:Kt + 6t) 



P(x) 



where 



where 



P(n\x) 



(27rA)" 



-1/2 



exp 



{n — nY 
2A 



X 



n 



A 



2xSt' 

K 



(29) 
(30) 

(31) 
(32) 



For the conditional state to be sharply peaked on a particular number n we require 
A < 1, that is 



8T6t > 1, 



(33) 



where F = ^ as above. If we now set St ~ Tn^^^/10 we arrive at the condition 

F > -f[N{nmax + 1) + (iV + l)nmax]. (34) 

This is the same as the fast-measurement condition arrived at previously. Therefore, 
phonon number quantum jumps between Fock states {|n) : n < n^ax} with lifetimes on 
the order of the typical lifetime r„ or greater can in principle always be resolved from the 
measurement record by integrating the homodyne photocurrent. This was illustrated 
above in figure [TJ 
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4-. 5. Conditional phonon number statistics 

It is also interesting to look at the conditional statistics of the Fock state distribution. 
We integrate equation ( |T9|) with k = 1007iV and / = 4007iV for a total time of 
30 / {'~fN) seconds, and bin the mean phonon number at each time step of the simulation 
to the nearest integer to calculate the distribution. Figure H] shows the results for bath 
temperatures of iV = 0.5 and N = 1. The deviations in the numerical data from thermal 
statistics may be due to the measurement bath preventing the mechanical mode from 
completely thermalising with the phonon bath. 
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Figure 4. Conditional statistics of the Fock state distribution when the system is well 
within the good-measurement limit, (a) N — 0.5, (b) N — 1. 



4-6. Experimental prospects 

We discuss here the experimental prospects for observing phonon number quantum 
jumps in relation to the membrane-in-the-middle setup employed in [H]. Assuming 
ground state cooling of the mechanical mode and a bath temperature much higher than 
the energy level spacing of the mode, the adiabatic and fast-measurement conditions 
respectively become k, and x^/^ ^ ^f"; where Q is the quality factor of the 

mechanical mode and T is the bath temperature. At T = 300mi^ the mechanical mode 
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has Q = 1.2 X 10^. Thus we find the requirements for observation of the jumps to be 
K ^ 3x lO'^s^^ and / k ^ 3x lO^s"^. The first condition is satisfied by the setup of |T3], 
which has k ~ 0.3 x lO^s^^. Satisfying the second condition, however, remains a great 
challenge. The largest coupling achieved in the setup of [14] is G ~ 3 x 10~^s~^. With an 
incident power of 10/iiy we find x ~ 10^s~^. This is far short of the fast-measurement 
requirement ^ 10^s~^, and significant further experimental effort will be required 
to make observation of phonon number quantum jumps feasible. Fabricating ultra-low 
absorption membranes may allow for higher incident powers, thereby enhancing a^. 
Reducing the mode frequency and the motional mass of the mode, and finding ways to 
further increase the frequency shift of the cavity can enhance G. Finally, lower bath 
temperatures and higher Q can reduce the required size of x- 

5. Conclusions 

Optomechanical systems with coupling that is quadratic in the mechanical displacement 
degree of freedom are promising candidates for monitoring phonon number quantum 
jumps in macroscopic mechanical oscillator modes. For homodyne detection on a 
strongly driven single-port cavity, we have established two conditions that must be 
satisfied in these systems in order to observe the quantum jumps; the adiabatic condition 
requires the damping rate k of the cavity to be much larger than the thermalization 
rate of the mechanical mode, and the fast-measurement condition requires the quantity 
X^ / to be much larger than the thermalization rate of the mechanical mode, where x 
is the effective optomechanical coupling strength. We have shown through numerical 
integration of the SME that satisfying these two constraints gives rise to phonon number 
quantum jumps in both the strong coupling ^ 1) and weak coupling < 1) 

regimes. 
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